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Abstract 

Motivated by measurements of the flavor singlet axial coupling constant of the nucleon in po- 
larized deep inelastic scattering we study the contribution of instantons to OZI violation in the 
axial-vector channel. We consider, in particular, the fi — oi meson splitting, the flavor singlet and 
triplet axial coupling of a constituent quark, and the axial coupling constant of the nucleon. We 
show that instantons provide a short distance contribution to OZI violating correlation functions 
which is repulsive in the /i meson channel and adds to the flavor singlet three-point function of 
a constituent quark. We also show that the sign of this contribution is determined by positivity 
arguments. We compute long distance contributions using numerical simulations of the instanton 
liquid. We find that the iso-vector axial coupling constant of a constituent quark is {Qj^q = 0.9 
and that of a nucleon is g\ = 1.28, in good agreement with experiment. The flavor singlet coupling 
of a quark is close to one, while that of a nucleon is suppressed, g\ = 0.77. However, this number 
is larger than the experimental value g\ = (0.28 — 0.41). 
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I. INTRODUCTION 



The current interest in the spin structure of the nucleon dates from the 1987 discovery by 
the European Muon Collaboration that only about 30% of the spin of the proton is carried 



by the spin of the quarks 3|- This results is surprising from the point of view of the naive 
quark model, and it implies a large amount of OZI (Okubo-Zweig-Iizuka rule) violation in 
the flavor singlet axial vector channel. The axial vector couplings of the nucleon are related 
to polarized quark densities by 

g\ = ^u- Ad, (1) 
g\ = Au + Ad-2As, (2) 
c/^ = AS = Am + Ad + As. (3) 

The first linear combination is the well known axial vector coupling measured in neutron 
beta decay, g\ = 1.267 ± 0.004. The hyperon decay constant is less well determined. A 
conservative estimate is g\ = 0.57 ± 0.06. Polarized deep inelastic scattering is sensitive 
to another linear combination of the polarized quark densities and provides a measurement 



of the flavor sing! 
(0.28-0.41), see 



et axial coupling constant g^. Typical results are in the range g\ 
21 for a recent review. 



Since g^ is related to the nucleon matrix element of the flavor singlet axial vector current 



many authors have specu 



ated that the small value of g^^ is in some way connected to the 



axial anomaly, see 

urn 

for reviews. The axial anomaly relation 

d'-Al = ^G;^G% + ^ 2m^g7z75g/ (4) 

implies that matrix elements of the flavor singlet axial current A° are related to matrix ele- 
ments of the topological charge density. The anomaly also implies that there is a mechanism 
for transferring polarization from quarks to gluons. In perturbation theory the nature of 
the anomalous contribution to the polarized quark distribution depends on the renormal- 
ization scheme. The first moment of the polarized quark density in the modified minimal 



subtraction (MS) scheme is related to the first moment in the Adler-Bardeen (AB) scheme 



by 



AE^ = AS^^ - Nf^^ACm, (5) 
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where AG is the polarized gluon density. Several authors have suggested that A'Eab is 
more naturally associated with the "constituent" quark spin contribution to the nucleon 
spin, and that the smallness of AS;^ is due to a cancellation between AT^ab and AC The 
disadvantage of this scheme is that AS^^ is not associated with a gauge invariant local 
operator |7]. 

Non-perturbatively the anomaly implies that = AS can be extracted from nucleon 
matrix elements of the topological charge density fi'^G^j^G^j,/ (327r^) and the pseudoscalar 
density mipi'-jr^ip . The nucleon matrix element of the topological charge density is not known, 
but the matrix element of the scalar density g^GI^^G^ is fixed by the trace anomaly. We 



have 



{N{p)\^,Gl^G%\N{p')) = Cs{q')mr,u{p)u{p'), (6) 

with Gs{0) = —1/b where b = 11 — 2Nf/3 is the first coefficient of the QCD beta function. 
Here, u{p) is a free nucleon spinor. Anselm suggested that in an instanton model of the QCD 
vacuum the gauge field is approximately self-dual, G^ = ±GG, and the nucleon coupling 
constants of the scalar and pseudoscalar gluon density are expected to be equal, Cs{0) ~ 
Cp(0) Q, see also Q- Using c/^ = NfGp{0) in the chiral limit we get ^ -Nf/b ^ -0.2, 
which is indeed quite small. 



A different suggestion was made by Narison, Shore, and Veneziano Narison et 

al. argued that the smallness of AS = is not related to the structure of the nucleon, 
but a consequence of the U{1)a anomaly and the structure of the QCD vacuum. Using 
certain assumptions about the nucleon-axial- vector current three-point function they derive 
a relation between the singlet and octet matrix elements, 

9'A = 9\^^fiU0)- (7) 

Here, = 93 MeV is the pion decay constant and Xtop(O) is the slope of the topological 
charge correlator 

Xtopiq') = J d'xe''^^ {Qtop{0)Qtop{x)) , (8) 

with Qtop{x) = (7^G^^G^^/(327r^). In QCD with massless fermions topological charge is 
screened and Xtap(O) = 0. The slope of the topological charge correlator is proportional to 
the screening length. In QCD we expect the inverse screening length to be related to the r]' 
mass. Since the r]' is heavy, the screening length is short and Xtopi^) is small. Equation ((Tj) 
relates the suppression of the flavor singlet axial charge to the large r]' mass in QCD. 



Both of these suggestions are very interesting, but the status of the underlying assump- 
tions is somewhat unclear. In this paper we would like to address the role of the anomaly 
in the nucleon spin problem, and the more general question of OZI violation in the flavor 
singlet axial-vector channel, by computing the axial charge of the nucleon and the axial- 
vector two-point function in the instanton model. There are several reasons why instantons 
are important in the spin problem. First of all, instantons provide an explicit, weak cou- 
pling, realization of the anomaly relation equ. Q and the phenomenon of topological charge 
screenin ej ll^ IfS^ . Second, instantons provide a successful phenomenology of OZI violation 
in QCD ^]. Instantons explain, in particular, why violations of the OZI rule in scalar me- 
son channels are so much bigger than OZI violation in vector meson channels. And finally, 
the instanton liquid model gives a very successful description of baryon correlation functions 



and the mass of the nucleon 
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This paper is organized as follows. In Sect. |n] we review the calculation of the anomalous 
contribution to the axial- vector current in the field of an instanton. In Sect. II I II and HVl we 
use this result in order to study OZI violation in the axial- vector correlation function and the 
axial coupling of a constituent quark. Our strategy is to compute the short distance behavior 
of the correlation functions in the single instanton approximation and to determine the large 
distance behavior using numerical simulations. In Sect.|3we present numerical calculations 
of the axial couplings of the nucleon and in Sect. IVII we discuss our conclusions. Some 
results regarding the spectral representation of nucleon three-point functions are collected 
in an appendix. 



II. AXIAL CHARGE VIOLATION IN THE FIELD OF AN INSTANTON 

We would like to start by showing explicitly how the axial anomaly is realized in the 
field of an instanton. This discussion will be useful for the calculation of the OZI violating 
part of the axial-vector correlation function and the axial charge of the nucleon. The flavor 
singlet axial-vector current in a gluon background is given by 

A^{x) = Tt[^,^,S{x,x)] (9) 

where S{x, y) is the full quark propagator in the background field. The expression on the 
right hand side of equ. Q is singular and needs to be defined more carefully. We will employ 
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a gauge invariant point-splitting regularization 



^T['y5l^lS{x,x)] = liniTr 



/ rx+e 



(10) 



In the following we will consider an (anti) instanton in singular gauge. The gauge potential 
of an instanton of size p and position 2; = is given by 



K-^%^nl. (u) 



Here, r]'^^ is the 't Hooft symbol and R characterizes the color orientation of the instanton. 
The fermion propagator in a general gauge potential can be written as 

^ A — m 

where x{x) is a normalized eigenvector of the Dirac operator with eigenvalue A, Ip'^\{x) = 
X'^xlx). We will consider the limit of small quark masses. Expanding equ. (fT^ in powers 
of m gives 

. _££(£MM + J2 + ,n J: MflM + o(,„^). (13) 

Here we have explicitly isolated the zero mode propagator. The zero mode \l/o was found by 
't Hooft and is given by 

Here, 0'^" = e"°/V2 is a constant spinor and 7± = (1 ±75)72 for an instanton/anti-instanton. 
The second term in equ. (fT^ is the non-zero mode part of the propagator in the limit m —>■ 



A^O ^ 

where Z)^ = — iA'^ and A±{x,y) is the propagator of a scalar field in the fundamental 
representation. Equ. (fT3j) can be verified by checking that 5*^^ satisfies the equation of 
motion and is orthogonal to the zero mode. The scalar propagator can be found explicitly 

A±(x, y) = Ao(x, y)^l= (l + ^'^^ " ^ ) -j^, (16) 

where Aq = l/(47r^A^) with A = x — y is the free scalar propagator. The explicit form of 
the non-zero mode propagator can be obtained by substituting equ. (fTHjl into equ. (|TH|l . We 



find 

^1 + g .^1 + ^ I V J 

X y \p + 

+ p2^^2 ^T • ■ ^^T ■ 7^± • yiT^ I (17) 

Here, So = — Z)(/(27r^A^) denotes the free quark propagator. As expected, tlie full non-zero 
mode propagator reduces to the free propagator at short distance. The linear mass term in 
equ. ()13|) can be written in terms of the non-zero mode propagator 

E ^^^m^ = / dhS^^ix,z)S^'iz,y) = -A±(x,y)7± - A^J{x,yh^, (18) 

where A±(a;,?/) is the scalar propagator and A^^(x,?/) = {x\{D'^ + a ■ Gl2)^^\y) is the 
propagator of a scalar particle with a chromo magnetic moment. We will not need the 
explicit form of A^ (x, y) in what follows. We are now in the position to compute the 
regularized axial current given in equ. ()10|) . We observe that neither the free propagator nor 
the zero mode part will contribute. Expanding the non-zero mode propagator and the path 
ordered exponential in powers of e we find 

Tr [1,1^ S{x, x)] = ± ,,^f';\,3 , (19) 

which shows that instantons act as sources and sinks for the flavor singlet axial current. We 
can now compare this result to the anomaly relation equ. (^. The divergence of equ. (|19p 
is given by 



^M-) = i^Sx^Si- (20) 



2p^ 4p^ - 2x 

^2 (-^2 _|_ p 

The topological charge density in the field of an (anti) instanton is 



We observe that the divergence of the axial current given in equ. ()19p does not agree with the 
topological charge density. The reason is that in the field of an instanton the second term 
in the anomaly relation, which is proportional to mip'y^ip, receives a zero mode contribution 
and is enhanced by a factor 1/m. In the field of an (anti) instanton we find 

2m^z75^ = T^^4t-^- (22) 



Taking into account both equ. ()2ip and ()22j) we find that the anomaly relation Q is indeed 
satisfied. 



III. OZI VIOLATION IN AXIAL- VECTOR TWO-POINT FUNCTIONS 

In this section we wish to study OZI violation in the axial-vector channel due to instan- 
tons. We consider the correlation functions 

W^,{x,y) = {j^{x)jM). (23) 

where is one of the currents 



where in the brackets we have indicated the mesons with the corresponding quantum num- 
bers. We will work in the chiral limit niu = frid 0. The iso- vector correlation functions 
only receive contributions from connected diagrams. The iso-vector vector (p) correlation 
function is 

inl),A^,y) = 2iPr),uix,y) = -2(Tr [^^S{x,yh,Siy,x)]). (25) 

The iso-singlet correlator receives additional, disconnected, contributions, see Fig. ^ The 
iso-singlet vector {u) correlator is given by 

iK),Ax,y) = 2(P-")^,(x,y) +4(P^-)^,(a;,y) (26) 

with 

{P^n,Ax,y) = (Tr [^,S{x,x)] Tr [i.S{y,y)]). (27) 

The axial-vector correlation functions are defined analogously. At very short distance the 
correlation functions are dominated by the free quark contribution Il\ = U^^^ = Uy = Uy ~ 
1/x^. Perturbative corrections to the connected correlators are 0{as{x)/7T), but perturbative 
corrections to the disconnected correlators are very small, 0{{as{x)/TT)^). In this section we 
will compute the instanton contribution to the correlation functions. At short distance, it is 
sufficient to consider a single instanton. For the connected correlation functions, this calcu- 
lation was first performed by Andrei and Gross see also jl^. Disconnected correlation 
function were first considered in {2^ and a more recent study can be found in 21 1. 



In order to make contact with our calculation of the vector and axial-vector three-point 
functions in the next section, we briefly review the calculation of Andrei and Gross, and 
then compute the disconnected contribution. Using the expansion in powers of the quark 
mass, equ. (fT^ . we can write 



(28) 



with 



/ T:>con\liU 
{-^V JO 



-Tr 
-2Tr 



(29) 
(30) 
(31) 



Using the explicit expression for the propagators given in the previous section we find 



and 



-/i,(x^A, + Xc,/\p)\ T 2e'^""^(/ij,A„?/^ - K/\px^)] (32) 



liihl {{x-y + p'^W' - {y^x' - x^y') T e^'^'^x^yp] (33) 



with K = l/(a;^ + p^), A = a; - ?/, and S''""'^ = gf^'g"!^ - gi'-'g^i^ + ^m/^^"^. Qur result agrees 
with [l^ up to a color factor of 3/2, first noticed in [22I, a '-' sign in front of the epsilon 
terms, which cancels after adding instantons and anti-instantons, and a '-' sign in front of 
the 2nd term in Bf^" . This sign is important in order to have a conserved current, but it does 
not affect the trace P^^^^. Summing up the contributions from instantons and anti-instantons 
we obtain 



(27r2)2A8 27T^ 



2S^ 



A4 



A^A" 



2S ■ A 



(S^A'' + A^S'' 



-S ■ Ag^"") + J^(A V - A^A'' - A^S'^ + A'^S'^) 



(34) 



This result has to be averaged over the position of the instanton. We find 



2a^'' = Y. I d^zA''^{x -z,y-z) 



± [X - z,y - z) = — 



Q2 



G{A\p) + 2G\A\p)g 



(9A2 



G{A\p) + 2G\A\p) 



(35) 
(36) 



with 



G'(A2,p) 



f loe; 



A 



1 



(37) 



9A2 A4 

and ^2 = AV(A2 + 4p2). The final resuh for the single instanton contribution to the 
connected part of the vector current correlation function is 



V ) - y 
where we defined 



2(P, 



con \ /^M 
V JO 



24 d 



7r2A2 aA2 



9A2 I A 



2 log 



24 rt^ 

-|^i^(A,0, (38) 



TT 



(39) 



The computation of the connected part of the axial-vector correlator is very similar. Using 
we observe that the only difference is the sign in front of B^'^ . We find 

1 



equs. 



6[P^'^)^^ = (^p^'^'ji^f^ _ 2(P^°")q'^ 



7r^ 



20A'G"' + 56G" 



(40) 



with G' given in equ. (j?7j) 

We now come to the disconnected part, see Fig. |21 In the vector channel the single in- 
stanton contribution to the disconnected correlator vanishes . In the axial- vector channel 
we can use the result for Tr[757^S'(x, x)] derived in the previous section. The correlation 
function is 



{prr{x,y) 



(41) 



Summing over instantons and anti-instantons and integrating over the center of the instanton 
gives 



'dis\av 
A ) 



27r2 dA^dA^ 



F{A,0 



and 



4p^ 



d 



dA^ 



+ A' 



dA^) ] 



F{A,0. 



(42) 



(43) 



We can now summarize the results in the vector singlet (u) and triplet (p), as well as 
axial- vector singlet (/i) and triplet (ai) channel. The result in the p and u channel is 



MM 



12 



r 24 
+ 2j dpnip)-^^FiA,0- 



TT^Af 



In the Oi, /i channel we have 
12 



vr^AS 
12 



TT 



4^6 



dpn{p) 
dpn{p) 



(20A2G" + 56G" 



20A^G" + 56G' 



16p^ 



2F' + A^F"^ 



(44) 

(45) 
• (46) 
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In order to obtain a numerical estimate of the instanton contribution we use a very simple 
model for the instanton size distribution, n{p) = no6{p — p), with p = 0.3 fm and uq = 
0.5 fm~^. The results are shown in Fig. |21 

We observe that the OZI rule violating difference between the singlet and triplet axial- 
vector correlation functions is very small and repulsive. We can also see this by studying 
the short distance behavior of the correlation functions. The non-singlet correlators satisfy 



(n^)'^'^ = (nf.)r |i + ^ / dp^ip) + •••}' (47) 

{^'Ar = (ni)r {l - ^'^'J dpnip) + ...}. (48) 

As explained by Dubovikov and Smilga, this result can be understood in terms of the 
contribution of the dimension d = 4 operators {g'^G'^) and {mqq) in the operator product 
expansion (OPE). The OZI violating contribution 

ij^T'r = {^\r' - = (^^) / dpn{p) (49) 

is of O(x^) and not singular at short distance. Our results show that it remains small 
and repulsive even if x > p. We also note that the sign of the OZI-violating term at 
short distance is model independent. The quark propagator in euclidean space satisfies the 
Weingarten relation 

^(x,?/)^ = 755'(?/,x)75. (50) 

This relation implies that Tr[S'(a;, x)7^75] is purely real. As a consequence we have 

lim |Tr[S'(a;, x)7^75]Tr[5'(?/, y)!^,!^]] > 0. (51) 

Since the path integral measure in euclidean space is positive this inequality translates into 
an inequality for the correlation functions. In our convention the trace of the free correlation 
function is negative, and equ. implies that the interaction is repulsive at short distance. 
The result is in agreement with the single instanton calculation. 

We can also study higher order corrections to the single instanton result. The two- 
instanton (anti-instanton) contributions is of the same form as the one-instanton result. An 
interesting contribution arises from instanton-anti-instanton pairs, see Fig. El This effect 

n 

was studied in .23] . It was shown that the instanton-anti-instanton contribution to the 

disconnected meson channels can be described in terms of an effective lagrangian 

2C 

C = —,{^l,l,i^Y (52) 
10 



with 

G= f dp,dp,{2^p,f{2^p,f'^^^ (53) 

where n{pi,p2) is the tunnehng rate for an instanton-anti-instanton pair and Tj^ is the 
matrix element of the Dirac between the two (approximate) zero modes. We note that this 
interaction is also repulsive, and that there is no contribution to the u channel. 

Numerical results for the vector meson correlation functions are shown in Figs. El and ^ 
unctions are obtained from Monte Carlo simulations of the instanton liquid 



The correlation 

as described in |24L |25|. We observe that OZI violation in the vector channel is extremely 
small, both in quenched and unquenched simulations. The OZI violating contribution to 
the /i channel is repulsive. In quenched simulations this contribution becomes sizable at 
large distance. Most likely this is due to mixing with an rj' ghost pole. We observe that the 
effect disappears in unquenched simulations. The pion contribution to the Oi correlator is 
of course present in both quenched and unquenched simulations. 

Experimentally we know that the p and u, as well as the Oi and /i meson, are indeed 
almost degenerate. Both iso-singlet states are slightly heavier than their iso- vector partners. 
To the best of our knowledge there has been only one attempt to measure OZI violating 
correlation functions in the vector and axial- vector channel on the lattice, see (26|. Isgur 
and Thacker concluded that OZI violation in both channels was too small to be reliably 
measurable in their simulation. 



IV. AXIAL VECTOR COUPLING OF A QUARK 

In this section we wish to study the iso-vector and iso-singlet axial coupling of a single 
quark. Our purpose is twofold. One reason is that the calculation of the axial-vector 
three-point function involving a single quark is much simpler than that of the nucleon, and 
that it is closely connected to the axial-vector two-point function studied in the previous 
section. The second, more important, reason is the success of the constituent quark model 
in describing many properties of the nucleon. It is clear that constituent quarks have an 
intrinsic structure, and that the axial decay constant of a constituent quark need not be 
close to one. Indeed, Weinberg argued that the axial coupling of a quark is {qaJq — 0-8 
Using this value of (5'i)Q together with the naive SU{6) wave function of the nucleon gives 
the nucleon axial coupling g\ = 0.8 • 5/3 ~ 1.3, which is a significant improvement over the 
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naive quark model result 5/3. It is interesting to study whether, in a similar fashion, the 
suppression of the flavor singlet axial charge takes place on the level of a constituent quark. 

In order to address this question we study three-point functions involving both singlet 
and triplet vector and axial-vector currents. The vector three-point function is 

{WyQQ)f{x,z,y) = {q"{x)V^-{z)q^iy)). (54) 

The axial-vector function (II^qq) is defined analogously. We should note that equ. (fH^ is 
not gauge invariant. We can define a gauge invariant correlation function by including a 
gauge string. The gauge string can be interpreted as the propagator of a heavy anti-quark, 
see Fig. IHl This implies that the gauge invariant quark axial-vector three-point function is 
related to light quark weak transitions in heavy-light mesons. 

The spectral representation of vector and axial- vector three-point functions is studied in 
some detail in the appendix. The main result is that in the limit that ^ ^ X4 the 
ratio 

Tr[(n^QQ) 37573] gA .ggx 
Tr[(n^QQ)474] ^ gv ^ ^ 

tends to the ratio of axial- vector and vector coupling constants, Qa/ov- We therefore define 
the following Dirac traces 

{Ut.QQrix,z,y) = Tr[(n»^QQ)V], (56) 
(W^Q^rix^z^y) = Tr[(n^QQ)'^757l- (57) 

As in the case of the two-point function the iso-triplet correlator only receives quark-line 
connected contributions, whereas the iso-singlet correlation function has a disconnected 
contribution, see Fig. El We find 

{Ul.QQr{x,z,y) = {P^"-Qr{x,z,y) = {Tt[S{x, z)rS{z,y)Y]) (58) 

and 

{U%Qr{x,z,y) = {P^"-Qr{x,z,y) - 2{P^^Qr{x, z,y) (59) 

with 

{P^QQr{x,z,y) = {TT[S{x,y)j^]TT[S{z,z)Y]), (60) 
as well as the analogous result for the axial- vector correlator. 
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In the following we compute the single instanton contribution to these correlation func- 
tions. We begin with the connected part. We again write the propagator in the field of 
the instanton as Szm + Sjqz + where Szm is the zero-mode term, S^z is the non-zero 
mode term, and 5*^ is the mass correction. In the three-point correlation function we get 
contribution of the type Sj^zSnz, SzuSm and SmSzM 

{PaJvQqY" = PnZNZ + CA/y {PzMm + PmZM) 

= TT[S''^{x,z)rS''^{z,y)Y] 

+CA/v^Tr[-vPo(^)*^W7'^(-A±(^,Z/)7±7l 
+CA/vm-A^{x,zh±m-Mz)^tiyW], (61) 

where the only difference between the vector and axial- vector case is the sign of 'ZMm' and 
'mZM' terms. We have ca/v = il for vector (axial vector) current insertions. The detailed 
evaluation of the traces is quite tedious and we relegate the results to the appendix iBl 

Our main goal is the calculation of the disconnected correlation function, which is related 
to OZI violation. In the single instanton approximation only the axial-vector correlator 
receives a non-zero disconnected contribution 

(PSgr = Tr[5(x,i/)757lTr[5(z,^)757l, (62) 

see Fig. We observe that the second trace is the axial-vector current in the field of an 
instanton, see equ. As for the first trace, it is easy to see that neither the zero mode part 
of the propagator nor the part of Snz proportional to the free propagator can contribute. 
A straight-forward computation gives 

Tr[^(x,y)757l = T 2 2 f . 4" ^ - 4" ^ • (63) 



Combined with equ. ()T9|) we obtain 

^ 2p^zVy^{x - yf 1 /^S^ _ 

y AQQ) 7r4(a; - yfx^y^z^ + p^f /(i + 4)(1 + 4) + + vV ' 



which has to be multiplied by a factor 2 in order to take into account both instantons and 
anti-instantons. 

Results for the vector and axial-vector three-point functions (n^QQ)'^^(r, r/2, 0) and 
(n^QQ)^'^(r, r/2, 0) are shown in Fig.|H| We observe that the vector and axial-vector correla- 
tion functions are very close to one another. We also note that the disconnected contribution 
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adds to the connected part of the axial- vector three-point function. This can be understood 
from the short distance behavior of the correlation function. The disconnected part of the 
gauge invariant three-point function satisfies 

\nn^{{U%Qr{x,z,y)} = lim {Tr[5(a;, a;)7375]Tr[5(^, ^)7375]} > 0. (65) 

This expression is exactly equal to the short distance term in the disconnected /i meson 
correlation function. The short distance behavior of the connected three-point function, on 
the other hand, is opposite in sign to the two-point function. This is related to the fact 
that the two-point function involves one propagator in the forward direction and one in 
the backward direction, whereas the three-point function involves two forward propagating 
quarks. A similar connection between the interaction in the /i meson channel and the flavor 



singlet coupling of a constituent quark was found in a Nambu-Jona-Lasinio model 
It was observed, in particular, that an attractive coupling in the /i channel is needed in 
order to suppress the flavor singlet ((^^ )q. 

The same general arguments apply to the short distance contribution from instanton- 
anti-instanton pairs. At long distance, on the other hand, we expect that lA pairs reduce 
the flavor singlet axial current correlation function. The idea can be understood from Fig. [71 
see Sy Q. In a. lA .....tion a left-handed valence up ,ua.l. e,nits a ngKt Kanded down 
quark which acts to shield its axial charge. We have studied this problem numerically, see 
Figs. ISHTUl We find that in quenched simulations the flavor singlet axial three-point function 
is significantly enhanced. This effect is analogous to what we observed in the /i channel and 
disappears in unquenched simulations. We have also studied the axial three-point function 
at zero three- momentum q = 0. This correlation function is directly related to the coupling 
constant, see App.|X| We find that the iso-vector coupling is smaller than one, {g\)Q — 0.9, 
in agreement with Weinberg's idea. The flavor singlet coupling, on the other hand, is close 
to one. We observe no suppression of the singlet charge of a constituent quark. We have 
also checked that this result remains unchanged in unquenched simulations. 

V. AXIAL STRUCTURE OF THE NUCLEON 

In this section we shall study the axial charge of the nucleon in the instanton model. We 
consider the same correlation functions as in the previous section, but with the quark field 
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replaced by a nucleon current. The vector three-point function is given by 

{^vNN)f{^.y) = {tm',{yW{^))- (66) 

Here, 7]°" is a current with the quantum numbers of the nucleon. Three-quark currents with 
the nucleon quantum numbers were introduced by loffe 3^. He showed that there are two 
independent currents with no derivatives and the minimum number of quark fields that have 
positive parity and spin 1/2. In the case of the proton, the two currents are 

?7i = eabc{u''C-ff,u^)j5-f^d'', r]2 = eabc{u°-Caf,yU^)-i^a^ud''. (67) 

It is sometimes useful to rewrite these currents in terms of scalar and pseudo-scalar diquark 
currents. We find 

r/i = 2 [,ai>c{u''Cd')^^u' - e,6,(M"C75f^')M^} , (68) 
m = ^ [eabc{u''Cd!')^^u' + eaUu''C^^d')u'} . (69) 

Instantons induce a strongly attractive interaction in the scalar diquark channel 
e'^'^'^^u^C'j^d'^) a consequence, the nucleon mainly couples to the scalar diquark 

component of the loffe currents 771^2- This phenomenon was also observed on the lattice js^. 
This result is suggestive of a model of the spin structure that is quite different from the 
naive quark model. In this picture the nucleon consists of a tightly bound scalar-isoscalar 
diquark, loosely coupled to the third quark js^. The quark-diquark model suggests that the 
spin and isospin of the nucleon are mostly carried by a single constituent quark, and that 

N Q 

9a - 9a- 

Nucleon correlation functions are defined by H^^(x) = (?7Q,(0)f/^(x)), where a,P are Dirac 
indices. The correlation function of the first loffe current is 



Ilafsix) = 2eatcea'i>'c' {[l.l^S^ (0, xhu76)^^ Tr [7^^r (0, xh^St' (0, x))' C\). (70) 

The vector and axial-vector three-point functions can be constructed in terms of vector and 
axial-vector insertions into the quark propagator, 

{rl)f{x,y) = Sf{0,y)^,Sf'{y,x), (71) 
{T^)f{x,y) = Sf{0,y)j,j,S'/{y,x). (72) 
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The three-point function is given by all possible substitutions of equ. (f7T|l and (f?^ into the 
two-point function. We have 

(^VNN)fix.y) = '^^abcea'b'c' 

+ 2<7^>(7,75^f(0,x)7.75)^^ TT[^,{Tl)f{x,yh^C{St''{0,x)rc] (73) 
7p75^r'(0, a:)7.75)^^ Tr [i,Sf{0, x)7.C(Sf (0, x)fC 



■ Tr 



9v^t.S'^''iy^y) + 4'yt^sf{y,y) ) 



where the first term is the vector insertion into the d quark propagator in the proton, the 
second term is the insertion into the uu diquark, and the third term is the disconnected 
contribution, see Fig. ^2 The vector charges of the quarks are denoted by gy. In the case 
of the iso- vector three-point function we have gy = 1, gy = —1 and in the iso-scalar case 

Vector and axial-vector three-point functions of the nucleon are shown in Figs. ^1 and 
IT^ In order to verify that the correlation functions are dominated by the nucleon pole 
contribution we have compared our results to the spectral representation discussed in the 
appendix, see Fig. El The nucleon coupling constant was determined from the nucleon 
two-point function. The figure shows that we can describe the three-point functions using 
the phenomenological values of the vector and axial- vector coupling constants. We have also 
checked that the ratio of axial- vector and vector current three-point functions is independent 
of the nucleon interpolating field for x > 1 fm. The only exception is a pure pseudo-scalar 
diquark current, which has essentially no overlap with the nucleon wave function. 

The main result is that the iso-vector axial-vector correlation function is larger than the 
vector correlator. The corresponding ratio is shown in Fig. El together with the ratio of 
q = correlation functions. We find that the iso-vector axial coupling constant is g\ = 
1.28, in good agreement with the experimental value. We also observe that the ratio of 
point-to-point correlation functions is larger than this value. As explained in the appendix, 
this shows that the axial radius of the nucleon is smaller than the vector radius. Taking 
into account only the connected part of the correlation function we find a singlet coupling 
g'^ = 0.79. The disconnected part is very small, g\{dis) = —(0.02 ± 0.02). Assuming that 
As ~ Au{dis) = Ad{dis) this implies that the OZI violating difference g\ — g^ is small. 
This result does not change in going from the quenched approximation to full QCD. 
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We have also studied the dependence of the results on the average instanton size, see 
Fig. We observe that there is a slight decrease in the iso-singlet coupling and a small 
increase in the iso-vector coupling as the instanton size is decreased. What is surprising 
is that the disconnected term changes sign between p = 0.3 fm and p = 0.35 fm. The 
small value g\{dis) = —(0.02 ± 0.02) obtained above is related to the fact that the phe- 
nomenological value of the instanton size is close to the value where g\{dis) changes sign. 
However, even for p as small as 0.2 fm the disconnected contribution to the axial coupling 
g\{dis) = —(0.05 ± 0.02) is smaller in magnitude than phenomenology requires. 

VI. CONCLUSIONS 

The main issue raised by the EMC measurement of the flavor singlet axial coupling is 
not so much why is much smaller than one - except for the naive quark model there is 
no particular reason to expect g\ to be close to one - but why the OZI violating observable 
9a ~ 9a is large. Motivated by this question we have studied the contribution of instantons 
to OZI violation in the axial-vector channel. We considered the /i — ai meson splitting, 
the flavor singlet and triplet axial coupling of a constituent quark, and the axial coupling 
constant of the nucleon. We found that instantons provide a short distance contribution 
which is repulsive in the /i meson channel and adds to the gauge invariant flavor singlet 
three-point function of a constituent quark. We showed that the sign of this term is fixed 
by positivity arguments. 

We computed the axial coupling constants of the constituent quark and the nucleon using 
numerical simulations of the instanton liquid. We find that the iso-vector axial coupling 
constant of a constituent quark is {g\)Q = 0.9 and that of a nucleon is g\ = 1.28, in good 
agreement with experiment. The result is also in qualitative agreement with the constituent 
quark model relation g\ = 5/3 ■ {gAjg- The flavor singlet coupling of quark is close to one, 
while that of a nucleon is suppressed, g^ = 0.77. However, this value is still significantly 
larger than the experimental value g^ = (0.28 — 0.41). In addition to that, we find very 
little OZI violation. As ~ Au{dis) ^ —0.01. We observed, however, that larger values of 
the disconnected contribution can be obtained if the average instanton size is smaller than 
the phenomenological value of p ^ 1/3 fm. 

There are many questions that remain to be addressed. In order to understand what 
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is missing in our calculation it would clearly be useful to perform a systematic study of 
OZI violation in the axial- vector channel on the lattice. The main question is whether the 
small value of is a property of the nucleon, or whether large OZI violation is also seen 
in other channels. A study of the connected contributions to the axial coupling constant 
in cooled as well as quenched quantum QCD configurations was performed in |36|]. These 
authors find g\{con) = Au{con) + Ad{con) ~ 0.6 in both cooled and full configurations. The 



disconnected term was computed by Dong et al. [37|. They find Au{dis) + Ad{dis) ~ —0.24. 

In the context of the instanton model it is important to study whether the results for 
obtained from the axial-vector current three-point function are consistent with calculations 
of g^ based on the matrix element of the topological charge density GG 38, 40, 41 1. 
It would also be useful to further clarify the connection of the instanton liquid model to 
soliton models of the nucleon 42|. In soliton models the spin of the nucleon is mainly due 
to the collective rotation of the pion cloud, and a small value for g'^ is natural j43|, |44|. 
The natural parameter that can be used in order to study whether this picture is applicable 
is the number of colors, Nc- Unfortunately, a direct calculations of nucleon properties for 
Nc > 3 would be quite involved. Finally it would be useful to study axial form factors of 
the nucleon. It would be interesting to see whether there is a significant difference between 
the iso-vector and iso-singlet axial radius of the nucleon. A similar study of the vector form 
factors was recently presented in 4^. 

Acknowledgments: We would like to thank E. Shuryak for useful discussions. This work 
was supported in part by US DOE grant DE-FG-88ER40388. 

APPENDIX A: SPECTRAL REPRESENTATION 
1. Nucleon Two-Point Function 

Consider the euclidean correlation function 

ufix) = {tm^x)), (Al) 

where ri°'{x) is a nucleon current and a is a Dirac index. We can write 

n^^(x) = ni(x)(x ■ 7)"^ + M^W^- (A2) 
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The functions 111^2(3;) have spectral representations 

roo 

Ui{x) = / dspi{s)D'{,/^,x), (A3) 
Jo 

roo 

U2ix) = / dsp2is)D{^/^,x), (A4) 
Jo 

where pi,2(s) are spectral functions and 

TTl 

D{m,x) = j^Ki{mx), (A5) 

2 

777 

D'{m,x) = - K^imx), (A6) 

are the euclidean coordinate space propagator of a scalar particle with mass m and its 
derivative with respect to x. The contribution to the spectral function arising from a nucleon 
pole is 

= |A^|(5(s - m\), p2{s) = \X%\mN5{s - m\), (A7) 

where A^r is the coupling of the nucleon to the current, {Q\r]\N{p)) — X^uip), and m^r is the 
mass of the nucleon. It is often useful to consider the point-to-plane correlation function 

K%^{t)^ Jd'xUfir,x). (A8) 

The integral over the transverse plane insures that all intermediate states have zero three- 
momentum. The nucleon pole contribution to the point-to-plane correlation function is 

K^ir) = ^(1 + 74)"''|AAr|'exp(-mivT). (A9) 



2. Scalar Three-Point Functions 



Next we consider three-point functions. Before we get to three-point functions of spinor 
and vector currents we consider a simpler case in which the spin structure is absent. We 
study the three-point function of two scalar fields and a scalar current j. We define 

U{x,y) = {mjiy)<Pix))- (AlO) 

The spectral representation of the three-point function is complicated and in the following 
we will concentrate on the contribution from the lowest pole in the two-point function of 
the field 0. We define the coupling of this state to the field and the current j as 

(O|0(O)|$(p)) = A, (All) 
($(p')|i(0)|$(p)) = F{q'), (A12) 
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where F{q^) with q = p—p' is the scalar form factor. The pole contribution to the three-point 
function is 

n(x, y) = X^ J d^zD{m, y + z)F{z)D{m, x-y-z), (A13) 

where D{m, x) is the scalar propagator and F{z) is the Fourier transform of the form factor. 
In order to study the momentum space form factor directly it is convenient to integrate over 
the location of the endpoint in the transverse plane and Fourier transform with respect to 
the midpoint 

d'x I Ae^'^(0(O)j(r/2,y)0(r,x)) = exp(-mr)F(g2). (AM) 

The correlation function directly provides the form factor for spacelike momenta. Maiani 
and Testa showed that there is no simple procedure to obtain the time-like form factor from 
euclidean correlation functions 4(|. 

Form factors are often parametrized in terms of monopole, dipole, or monopole-dipole 
functions 



FM{q') = Fm{0)-/^, (A15) 



Fniq') = Fn{0) ( j,^^ ) , (A16) 

i2 



Q"^ + m 



V , 



W) = fMo(0)^(^^J . (A17) 

with = —q^. For these parametrization the Fourier transform to euclidean coordinate 
space can be performed analytically. We find 

Fm{x) = mlD{x,mv) (A18) 
Fd{x) = ml. l^~D\x, my) - D{x, my)) (A19) 

Fmd{x) = o ^ % ^ ^ :^{D{x,mi) - D{x,m2)) 

m2 — mf I — mf 



+ ^(^D'{x,m2) + D{x,m2))\ . (A20) 



We also consider three-point functions involving a vector current j^. The matrix element is 

mp')\j,mHp)) = q,F{q')- (A21) 

The pole contribution to the vector current three-point function is 

n^(x, y) = >^^ I d^zD{m, y + z)z^F'{z)D{m, x-y-z), (A22) 
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with F'{z) — dF{z)/dz and = ^t^/\^\- For the parameterizations given above the deriva- 
tive of the coordinate space form factor can be computed analytically. We get 

F'^ix) = m^yD'{x,mv) (A23) 
F'd(^) = -'^D{x,mv) (A24) 

Fmd{x) = -^^{-^^{D'{x,m,)-D'{x,m,)) 

+ |D(a;,m2)|. (A25) 



3. Three-Point Functions involving nucleons and vector or axial-vector currents 

Next we consider three-point functions of the nucleon involving vector and axial- vector 
currents. The vector three-point function is 

{^VNN)fM = mO)V^^{y)f{x)). (A26) 

The axial-vector three-point function is defined analogously. The nucleon pole contribution 

involves the nucleon coupling to the current 1] and the nucleon matrix element of the vector 
and axial vector currents. The vector current matrix element is 



Fiiq'h, + ^F,iq')a,^q'' 



\<P). (A27) 



{N{p')\V^^\N{p))=u{p') 
where the form factors Fi 2 are related to the electric and magnetic form factors via 



GE{q') = F^{q') + ^^F,{q'), (A28) 
GM{q') = F^{q') + F2{q'). (A29) 

The axial-vector current matrix element is 



{N{p')\A;\N{p))=u{p') GA{q'h, + 7^Gp{q'){p'-p\ 



l,-u{p) (A30) 



2M 

where Ga,p are the axial and induced pseudo-scalar form factors. 

We are interested in extracting the vector and axial- vector couphng constants gv — -^i(O) 
and qa — Ga(0). In order to determine the vector coupling gy we study the three-point 

function involving the four-component of the vector current in the euclidean time direction. 
For simplicity we take y = x/2. We find that 

iUy^Nr/{x, x/2) = + n^yiVAr(T)(74)"^ (A31) 
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where = i^^^)- The two independent structures Hy^jy are given by 



^ ^ '^^^mD'{xi)D{x2) + ^-^mD{xi)D'{x2) 



2xi 

X1X2 2m 
- Ayj + Af 



2xi 



(A32) 



4X1^2 



D\xi)D\x2) + m2D(xi)D(a;2) 



Fi{y) 



-D\x,)D{x2) + -D'(x2)L'(xi) 

Xl X2 



EM. 

2m 



(A33) 



where -^1,2(2/) are the Fourier transforms of the the Dirac form factors Fi_2(g^), and we have 
defined Xi = (y, t/2 + y^) and ^2 = {-y, t/2 - y^). 

In order to extract the axial-vector coupUng we study three-point functions involving 
spatial components of the axial- vector current. We choose the three-component of the 
current and again take y = x/2 with x^ = (0, r). We find 



{UANN)f{x, x/2) = nW(r)(75)'^^ + nW(T)(7375)"'' + nW(r)(737475)'^^ (A34) 



with 



1 



D {xi)D{x2) - 



n 



ANN 



V) 



d'y 



2xi 

t2 + 8y| - 4y2 



2X2 



D{xi)D'{x2) 



Ax 1X2 



D\xi)D'{x2) +m^D{xi)D{x2) 



+ 



\y\ 



771 Til 

-D'{xi)D{x2) + -D'{x2)D{xi] 

Xi X2 



G'pjy) 
2m 



Ga(2/),(A35) 
GA{y) 

(A36) 



where GA,p{y) are the Fourier transforms of the nucleon axial and induced pseudo-scalar 
form factors. 

These results are quite complicated. The situation simphfies if we consider three-point 
functions in which we integrate all points over their location in the transverse plane. The 
vector three-point function is 

j d'x I d'y{UvNN)2''{r,x;T/2,y) = ^(1 + 74)"'^|A^|2 exp(-m^T), ^^37) 

where gv — -^i(O) is the vector coupling. Note that the three-point function of the spatial 
components of the current vanishes when integrated over the transverse plane. The axial- 
vector three-point function is 

J d'x J d'y(UANN)f(r,x;r/2,y} = ^ ((1 + ^4)^3^5)"/^ |A^|2 exp(-mArT), (A38) 
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where qa = G'^(O) is the axial- vector couphng. In the case of the axial- vector current the 
three-point function of the time component of the current vanishes when integrated over 
the transverse plane. This is why we consider three-point function involving the spatial 
components of the axial current. 

4. Phenomenology 

In Fig.^lwe show the nucleon pole contribution to the vector and axial- vector three-point 
functions. We have used the phenomenological values of the isovector coupling constants 

G^(0) = 1, Gm(0) = 4.7, 

G^(0) = 1.25, Gp(0) = ^gA. 

We have parametrized Ge,v and Ga by dipole functions with my = 0.88 GeV and = 1.1 
GeV. The induced pseudoscalar form factor is parametrized as a pion propagator multiplied 
by a dipole form factor with dipole mass rriA- 

We observe that at distances that are accessible in lattice or instanton simulations, x ~ 
(1 — 2) fm, the typical momentum transfer is not small and the correlation function is 
substantially reduced as compared to the result for a point-like nucleon. We also observe 
that the F2 and Gp form factors make substantial contributions. Fig. [TBI shows that the 
ratio of the axial- vector and vector correlation functions is nevertheless close to the value for 
a point-like nucleon, qa/ Qv — 1-25. We observe that the ratio of point-to-point correlation 
functions approaches this value from above. This is related to the fact that the axial radius 
of the nucleon is smaller than the vector radius. As a consequence, the point-to-point 
correlation function at finite separation r "sees" a larger fraction of the axial charge as 
compared to the vector charge. 

APPENDIX B: INSTANTON CONTRIBUTION TO QUARK THREE-POINT 
FUNCTIONS 

In this appendix we provide the results for the traces that appear in the single instanton 
contribution to the quark three-point function. Our starting point is the expression 

(Bl) 
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see equ. fl61|) . Due to the Dirac structure of the non-zero mode part of the propagator, the 
NZNZ term is the same for both the vector and axial-vector correlation functions. It has 
4 parts stemming from combinations of the two terms of non-zero mode propagator equ. fTTI 



^NZNZii 



pfiu 

^NZNZi2 



Ti\x-z)\z-yY 
H{x, z, y){x — z)"'y^ 



1 + ^ 



p I X ■ z z ■ y 



+ 



\ x^ 



p^x ■ y 

x^z'^y^ 



(B2) 



27r4(x - z)\z - yfz'^y'^ 
p\z - yY^' 



x^ 



X 



^NZNZ2i 



p2 + z'^ 
H{x,z,y){z - vYxf^ 



p^{z - yY 



p2 + y2 



Te- 



rn 



[I'^ao _|_ i_,,ao^ c/3o-o-oao i r gp 



p2(x - zY° 



^NZNZ22 



p^ + 

2;, y)x°^°y^ 



p2 + z"^ 



(B4) 



47r^(x — — yYz'^y'^x'^ (p^ + 2;^) 



where 



(p2 + y2 



H{x,z,y) 



(p^ -|- 



1 + 



-1 



1 + 



x^ 



(B5) 



(B6) 



and the Dirac trace T± is defined as 



r^[p,2,3,4,5,6] = [(^A^2^3456 _^;.3^2456 ^^M^2356 _^m5^2346 ^^;.6^: 



2345 



where 2,3,... is short for p2, ^3, • • •• The SzAiSm term is easily seen to be 



(B7) 



pfiu 
^ZMm 



^^{z~yY ./(l + fi)(l + $) 



T:f[ao,o-i,a,/?o,p, i^] 



^(Ti,CT _|_ P ^ y ^^aiualS _j_ ^0-1 era/? ^ 



(B8) 



with y9(x) = pI{tx^{x^ + p^Y^"^)- The SmSzM term is obtained similarly 
^rzM = Tr[(-A±(x,^)7±)7'^(-^oW<(2/))7l 



X 



P2 
^ / ga-io-a/3 _|_ 

^,2^2 



(B9) 
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For all the above formulas, we need to add instanton and anti-instanton contribution and 
integrate over the position of the instanton, which was suppressed above. As usual, for an 
instanton at position zi we have to shift z according to x ^ {x — zj), etc. 
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FIG. 1: Quark line diagrams that contribute to the vector and axial- vector two-point function 
in the iso-vector and iso-singlet channel. The solid lines denote quark propagators in a gluonic 
background field. The two diagrams show the connected and disconnected contribution. 




FIG. 2: Instanton contributions to the disconnected axial- vector correlation function. The 
left panel shows the single-instanton (non-zero mode) contribution. The right panel shows the 
instanton-anti-instanton (fermion zero mode) contribution. 
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FIG. 3: Correlation functions in the p,uj,ai and /i channel. All correlation functions are nor- 
malized to free field behavior. The data points show results from a numerical simulation of the 
random instanton liquid. The dashed lines show the single instanton approximation. 
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FIG. 4: Correlation functions in the p,uj,ai and /i channel. All correlation functions are nor- 
malized to free field behavior. The data points show results from unquenched simulations of the 
instanton liquid model. 
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FIG. 5: Quark line diagrams that contribute to the vector and axial-vector three-point function of 
a constituent quark. The solid lines denote quark propagators in a gluonic background field. The 
two diagrams show the connected and disconnected contribution. 




FIG. 6: Physical interpretation of the gauge invariant axial- vector three-point function of a quark 
in terms of a weak light-quark transition in a heavy-light Qq meson. 
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FIG. 7: Instanton contributions to the disconnected axial-vector three-point correlation function 
of a quark. The left panel shows the single-instanton (non-zero mode) contribution. The right 
panel shows the instanton-anti-instanton (fermion zero mode) contribution. 
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FIG. 8: Axial and vector three-point functions of a quark as a function of the separation between 
the two quark sources. The correlation functions are normalized to free field behavior. The data 
points show results from numerical simulations of the instanton liquid and the dashed lines show 
the single instanton approximation. 
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FIG. 9: Axial and vector three-point functions of a quark as a function of the separation between 
the two quark sources. The data points show results from an unquenched instanton simulation. 
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FIG. 10: Ratio of axial- vector to vector correlation functions of a constituent quark calculated 
in the instanton liquid model. The open points show point-to-point correlation functions while 
the solid point is the zero momentum (point-to-plane) limit. The figure shows the iso-vector and 
iso-singlet correlation functions. 
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FIG. 11: Quark line diagrams that contribute to the axial- vector three-point function of the proton. 
The solid lines denote quark propagators in a gluonic background field. The lines are connected 
in the same way that the Dirac indices of the propagators are contracted. The iso-vector and 
iso-singlet correlation functions correspond to g\ = —gj^ = 1 and = gj^ = 1, respectively. The 
disconnected diagram only contributes to the iso-scalar three-point function. 



32 




0.5 1 1.5 

X [fm] 



FIG. 12: Vector, axial-vector three-point functions of the nucleon and nucleon two-point function 
calculated in the instanton liquid model. All correlation functions are normalized to free field 
behavior. The results are compared to a simple pole fit of the type discussed in the appendix. 
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FIG. 13: Ratio of axial-vector to vector correlation functions of the nucleon calculated in the 
instanton liquid model. The open points show point-to-point correlation functions while the solid 
point is the zero momentum (point-to-plane) limit. The figure shows the iso-vector, connected 
iso-singlet, and full iso-singlet axial-vector correlation functions. 




FIG. 14: Axial coupling constants of the nucleon as a function of the instanton size p with the 
instanton density fixed at {N/V) = lfm~^. We show the iso-vector, connected iso-singlet, and full 
iso-singlet axial coupling constant. 
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FIG. 15: Nucleoli pole contribution to the vector (upper panel) and axial-vector (lower panel) 
nucleon three-point function. The solid line shows the complete results, the dashed line is the 
contribution from the Fi and Ga form factors only, and the dash-dotted line corresponds to a point- 
like nucleon. We have used a nucleon coupling constant A = 2.2 fm~'^ as well as phenomenological 
values for the form factors and coupling constants. 
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FIG. 16: Ratio of the phenomenological parameterizations of the axial-vector and vector three- 
point functions. We have added a short distance continuum contribution to the nucleon pole terms. 
The curves are labeled as in the previous figure. Note that both the solid and the dashed line will 
approach = 1-25 as x — > oo. Also note that the solid line is in very good agreement with the 
instanton calculation shown in Fig. 1131 
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